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The sand fly Lutzomyia longipalpis is the most important vector of American Visceral 
Leishmaniasis. Adults are phytophagous (males and females) or blood feeders (females 
only), and larvae feed on solid detritus. Digestion in sand fly larvae has scarcely been 
studied, but some glycosidase activities putatively involved in microorganism digestion 
were already described. Nevertheless, the molecular nature of these enzymes, as the 
corresponding genes and transcripts, were not explored yet. Catabolism of microbial 
carbohydrates in insects generally involves p-1 ,3-glucanases, chitinases, and digestive 
lysozymes. In this work, the transcripts of digestive p-1 ,3-glucanase and chitinases were 
identified in the L. longipalpis larvae throughout analysis of sequences and expression 
patterns of glycoside hydrolases families 1 6, 1 8, and 22. The activity of one i-type lysozyme 
was also registered. Interestingly, this lysozyme seems to play a role in immunity, rather 
than digestion. This is the first attempt to identify the molecular nature of sand fly larval 
digestive enzymes. 
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INTRODUCTION 

Phlebotomines are psychodid dipterans distributed over almost 
all faunal regions of the world and particularly abundant in 
temperate and inter-tropical regions. There are more than 
600 described species to date. However, only 30 of these 
are epidemiologically important for pathogen transmission. 
Phlebotomine sand flies are able to transmit several diseases such 
as Leishmaniasis, which occur in approximately 98 countries and 
affects mOlions of people each year (WHO, 2014), as well as 
Bartonellosis (Tsai et al., 2011) and arbovirosis (Amaro et al., 
2012). 

Leishmaniasis are zoonosis caused by infection with 
Leishmania genus parasites and these are divided clinically 
in three forms (cutaneous, mucocutaneous, and visceral) 
depending on the parasite species. In the New World, Visceral 
Leishmaniasis is a severe systemic disease caused by Leishmania 
infantum (syn. L. chagasi) and transmitted by phlebotomines of 
the species Lutzomyia longipalpis (Romero and Boelaert, 2010; 
HarhayetaL, 2011). 

Like other Diptera, phlebotomines are holometabolous insects 
with adult and larval stages which exploit different food sources. 
Male and female adults feed on plant sap (or blood, in the case of 
females) (Brazil and Brazil, 2003; Scares and Turco, 2003) while 
larvae grow on decaying organic matter, mainly of plant origin. 



or animal feces. Wermelinger and Zanuncio (2001) found that 
L. longipalpis and L. intermedia larvae improved development 
until pupation when fed a variety of the diet described by Young 
et al. (1981), constituted by humus containing an intense pro- 
liferation of fungi. The authors suggested that this composition 
mimicked the larval substrate in the wilderness. However, despite 
these observations in the laboratory, little is known about the 
natural habitat of phlebotomine larvae (Feliciangeli, 2004). 

Differences between larval and adult feeding habits of sand 
flies correlate with anatomical, physiological, and biochemical 
variations within the phlebotomine digestive tract. Therefore, it is 
not accurate to generalize about digestion throughout the devel- 
opment of these insects. Although some authors have shown the 
presence of both proteases and carbohydrases in sand flies (Dillon 
and El Kordy, 1997; Gontijo et al., 1998; Jacobson and Schlein, 
2001; Ramalho-Ortigao and Traub-Csekd, 2003; Do Vale et al., 
2007; Telleria et al, 2007, 2010; Sant'Anna et al, 2009), the role 
of digestible carbohydrases in phlebotomines and their potential 
participation in the digestion of microorganisms has not been 
systematically studied. 

Moraes et al. (2012) described several glycosidases in larvae of 
L. longipalpis which may be related to the digestion of bacteria 
and fungi. This suggested that sand fly larvae have detritivo- 
rous feeding habits. The authors also assayed enzymatic activity 
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of P-glucanases, chitinases, and lysozymes and confirmed that 
P-l,3-glucanases exhibited the highest specific activities in the 
luminal contents of the larval gut. 

P-l,3-glucanases (PGlu) are glycoside hydrolases (GH) 
responsible for the hydrolysis of f5-l,3-glucans which are found in 
the cell walls of fungi (Bartnicki-Garcia, 1968; Gorin and Spencer, 
1968; Sietsma and Wessels, 1981) and phloem of higher plants as 
caUose (Bade et al., 1998). These enzymes are commonly found in 
the gut and salivary glands of insects such as Spodopterafrugiperda 
(Bragatto et al., 2010), Tenebrio molitor (Genta et al., 2009), and 
Periplaneta americana (Genta et al, 2003). 

Besides P-l,3-glucanases, insects also express fi-glucan recog- 
nition proteins (PGRPs), P-glucan binding proteins (GBP) and 
gram-negative bacteria binding proteins (GNBP) responsible 
for triggering the innate immune response through recogni- 
tion of pathogen-associated molecular proteins (PAMPs) such 
as lipopolysaccharide (LPS), peptidoglycan, and P-l,3-glucan 
(Royet, 2004). These proteins are found mainly in hemolymph 
and fat body of insects (Ma and Kanost, 2000; Wang et al, 2005; 
Sun et al., 2011). Both 3-glucanases and P-glucan-binding 
proteins found in insects have been structurally assigned to fam- 
ily 16 of glycoside hydrolases (GH16) (Genta et al, 2009; Bragatto 
etal, 2010). 

Chitinases (Chit) are enzymes able to hydrolyze chitin, a lin- 
ear polymer of N-acetylglucosamine (GlcNAc) present in the cell 
wall of fungi as well as in insect structures such as cuticle, trachea, 
peritrophic matrix, and gut (Arakane and Muthukrishnan, 2010). 
All insect chitinases belong to family 18 of glycoside hydrolases 
(GH18) and are responsible for cuticle turnover, as well as diges- 
tion and degradation of the peritrophic matrix during molting 
(Zhang et al, 2011). GH18 are encoded by several groups of genes 
which are constituted by a multi-domain structural organization 
that includes: (1) 1-5 catalytic domains; (2) 0-7 cysteine-rich 
chitin-binding domains (CBD), and (3) serine/threonine-rich 
linker regions that can be heavily glycosylated (Merzendorfer and 
Zimoch, 2003; Arakane and Muthukrishnan, 2010). 

Chitinases have been studied in different insect orders, in 
species such as Anopheles gamhiae (Shen and Jacobs-Lorena, 1997; 
Zhang et al, 2011), Aedes aegypti (De la Vega et al., 1998), 
Manduca sexta (Kramer et al., 1993), Bombyx mori (Kim et al., 
1998; Zhang et al, 2011), Hyphantria cunea (Kim et al., 1998), 
Drosophila melanogaster (Zhu et al., 2008), Ostrinia nubilalis 
(Khajuria et al, 2010), Tribolium castaneum (Zhu et al., 2008), 
Tenebrio molitor (Royer et al, 2002; Genta et al, 2006), and 
Phlebotomus papatasi (Ramalho-Ortigao et al., 2005; Coutinho- 
Abreu et al., 2010). Ramalho-Ortigao and Traub-Cseko (2003) 
isolated and characterized cDNA encoding a chitinase from the 
gut tissue of adult female L. longipalpis (Llchitl), which seems 
to be involved in the degradation of the peritrophic membranes 
and Leishmania migration toward the midgut epithelium during 
blood digestion. 

Lysozymes (Lys) are glycoside hydrolases belonging to the gly- 
coside hydrolases family 22 (GH22), whose function is to catalyze 
hydrolysis of glycosidic bonds between N-acetylmuramic acid 
(NAM) and N-acetylglucosamine (NAG). These molecules are 
constitutive components of the peptidoglycan layer of bacterial 
cell walls (JoUes and Jolles, 1984). Lysozymes are widely found in 



various organisms and are divided into six major groups. C-type 
lysozymes are probably the best studied and are found in several 
vertebrate and invertebrate taxa, including the vast class Insecta 
(Bachali et al., 2002; Harikrishnan et al, 2011). In addition to 
c-type lysozymes, other lysozymes have been studied, such as i- 
type lysozymes. This new class of lysozymes shows similarity with 
destabilases found in annelids (Zavalova et al., 2000), molluscs 
(Ito et al, 1999), and insects (Paskewitz et al., 2008). Zavalova 
et al. (2000) showed activity of i-type lysozyme in the cell wall of 
Micrococcus lysodeikticus and inhibition of such activity by anti- 
destabilase. However, Paskewitz et al. (2008) found no activity of 
i-type lysozyme in Anopheles gambiae. 

An important hindrance for the study of sand fly enzymes 
is the diminute size of these animals, and their laborious and 
time-consuming maintenance of colonies, which hardly affords 
samples with amounts of protein enough for activity screen- 
ings. For this reason, the traditional strategies of enzyme purifi- 
cation, characterization, and sequencing, which were used in 
the study of other insect digestive glucanases and chitinases 
(Genta et al., 2009), are not applicable to this insect model. 
Anterior work from our group (Moraes et al, 2012) showed 
the presence of fi-l,3-glucanase, chitinase, and lysozyme activ- 
ities in the gut of sand fly larvae, but there are no description 
of any protein or DNA sequence related to digestive enzymes 
in these insects. In this work, we studied the expression pat- 
tern of different genes from families GH16, GH18, and GH22 
in L. longipalpis and correlated two of them to the digestion 
of larvae. To our knowledge, this is the first report on the 
molecular identification of digestive enzymes in phlebotomine 
larvae. 

MATERIALS AND METHODS 

SAND FLY REARING AND LARVAE FEEDING 

Insects used in all experiments were insectary-reared L. longi- 
palpis from a colony originally started from individuals 
from Jacobina, Brazil, and maintained at the Laboratory of 
Insect Biochemistry and Physiology (Oswaldo Cruz Institute, 
FIOCRUZ). Adult sandflies were fed a 70% sucrose solution (w/v) 
ad libitum. Adult females were blood-fed with anesthetized ham- 
sters (ketamine, 200mg/kg) to trigger egg development. After 
oviposition, eggs were collected and reared to preserve the colony. 
All larval instars were fed a crushed mixture of rabbit feces, rabbit 
chow, and garden soil. Third and fourth instars were supple- 
mented with a mbcture of white soy protein (bran) and cereal 
flakes (Neston) (1:1). 

IN SILICO MINING OF L. LONGIPALPIS ESTs LIBRARY FOR GLYCOSIDE 
HYDROLASES 

To identify sequences of families 16 (P-1,3-Glucanases, GBP), 
18 (chitinases), and 22 (lysozymes), glycoside hydrolases from 
different insect orders were identified and retrieved from CAZy 
(Carbohydrate-Active Enzyme database, http://www.cazy.org/). 
Subsequently, GH sequences obtained from CAZy were employed 
to perform a TBLASTX search (http://blast.ncbi.nlm.nih.gov/) to 
find similar sequences in a Lu. longipalpis EST library (Dillon 
et al., 2006) at the Sanger Institute website (http://www.genedb. 
org/Page/parasite Vectors). 



Frontiers in Physiology | Irntegrative Physiology August 2014 | Volume 5 | Article 276 | 2 



Moraes et al. 



Sand fly larvae glycoside hydrolases 



ANALYSIS OF GH SEQUENCES FROM THE L. LONGIPALPIS EST LIBRARY 

Sequences retrieved from the L. longipalpis EST library were 
translated using the Translate tool (http://web.expasy.org/ 
translate/) and compared against the non-redundant pro- 
tein database from the National Center for Biotechnology 
Information (NCBI). Analysis ofLw. longipalpis ESTs similarity to 
GH families 16, 18, and 22 as well as prediction of ORE integrity 
was performed using the BLASTP tool. Translated sequences were 
further analyzed to identify signal peptides, O-type glycosylation, 
N-type glycosylation, and functional domains using PeptidelP 
Server 4.0 (Petersen et al., 2011; http://www.cbs.dtu.dk/services/ 
SignalP/), NetOGlyc Server 3.1 (lulenius et al, 2005; http://www. 
cbs.dtu.dk/services/NetOGlyc/), NetNGlyc Server 1.0 (http:// 
www.cbs.dtu.dk/services/NetNGlyc/), and database PEAM 26.0 
(Punta et al, 2012; http://pfam.sanger.ac.uk/), respectively. 
Alignments were performed using ClustalW (http://ebi.ac.uk/ 
Tools/msa/clustalw2/) (Hall, 1999). 

PHYLOGENETIC ANALYSIS OF L LONGIPALPIS GH16, GH18, AND GH22 

Neighbor-joining phylogenetic trees for each gene family were 
constructed using insect protein sequences with PEAM domains 
GH16_beta_GRP (CD02179, P-l,3-glucanases and P-glucan 
binding proteins), Glyco_hydro_18 (PE00704, chitinases), and 
LYS (PE00062) and Destabilase (CD05497) for lysozymes. Trees 
were generated using MEGA5.05 (Tamura et al., 2011). Bootstrap 
values were set at 5000 replications. 

DISSECTIONS 

4th instar larvae of Lu. longipalpis were rinsed in 200 [xL of ster- 
ile 0.15 M NaCl solutions and anesthetized on ice. Larvae were 
dissected and the following structures removed: head (includ- 
ing foregut), midgut and hindgut (M-l-H) and carcass (rest of 
the body). Matching tissues were pooled in groups of 5 and 
then transferred to polypropylene vials containing 50 (iL of TRI 
Reagent® (Sigma). Samples were flash-frozen and kept at — 80°C 
until further RNA extraction. 

RNA EXTRACTION AND cDNA SYNTHESIS 

RNA was extracted from entire insects, head, M-l-H and car- 
cass tissue samples, according to Diaz-Albiter et al. (2011). 
After extraction, total RNA was quantified using Nanodrop® 
(NanoDrop Technologies, Wilmington, USA). RNA was reverse- 
transcribed to cDNA using Superscript III Eirst-Strand Synthesis 
System (Invitrogen, San Diego, CA) following the manufacturer's 
protocol. cDNA was quantitated using Nanodrop and normalized 
to a concentration of 50 ng/|xl. 

TISSUE-SPECIFIC EXPRESSION OF pGLU, GBP, Chit, AND Lys 

Polymerase Chain Reaction (PGR) and multiplex PGR were used 
to assess tissue-specific expression of fi-Glu and GBP, Ghit, and 
Lys in fourth instar larvae. Eor multiplex PGR, three specific 
primers were included in the same reaction, which allowed simul- 
taneous amplification of more than one gene. Primer combina- 
tions were as follows: (1) llb04, LamS2, 24g06, and Ribo60; (2) 
96h07 and 154bl2; and (3) 123b01 and 18f06. Eor sequences 
88dl2 and 14b06, PGRs were performed using only a pair of 
primers. AH primers have listed in Table SI in Supplementary 



Material. Amplification reactions were performed in a total vol- 
ume of 20 |xL containing 50 ng cDNA, 0.5 (iM of each primer, 
1 X PGR reaction buffer, 0.2 mM each dNTP, 1.5 mM magnesium 
chloride and 0.025 U GoTaq® DNA polymerase (Promega). The 
parameters for PGR were: incubation at 94° G for 2 min followed 
by varying cycles of 94° G for 15 s, 55° G for 30 s, 72° G for 1 min 
and a final incubation (extension) of 72° G for 5 min. The num- 
ber of cycles vary depending on the combination of primers used 
in the PGR and PGR-Multiplex reactions, as shown in Table S2 
in the Supplementary Material. PGR products were analyzed by 
gel electrophoresis using agarose 2% (w/v), stained with ethidium 
bromide (0.5 |xg/mL). Expression patterns were obtained by mea- 
suring the band intensity by densitometry using ImageJ software 
and then calculating relative expressions against a constitutive 
gene (AM088777, 60S ribosomal protein L3). 

STATISTICAL ANALYSIS 

All statistical analyses were performed using GraphPad Prism 5.0 
for Windows, (San Diego Galifornia USA). Student's f-test for 
paired data was used. Comparison among and between groups 
was assessed using One- Way analysis of variance (ANOVA) with 
a post-hoc Tukey multiple comparison test. Results are expressed 
as the group mean ± s.e.m. Significance was considered when 
p < 0.05. 

RESULTS 

IN SILICO ANALYSIS OF GHs SEQUENCES RETRIEVED FROM THE 
L. LONGIPALPIS EST LIBRARY 

A total of 206 sequences of glycoside hydrolases families 16, 
18, and 22 from different insect species were found in GAZy, 
85 of which belonged to family GH16, 47 to GH18, and 
74 to GH22. They were employed to perform an in silico 
search targeting these three protein families in an L. longi- 
palpis ESTs library (http://www.genedb.org/Page/parasiteVectors, 
Dillon et al., 2006). A total of 9 sequences were retrieved, 3 
corresponding to clones of GH16 (identifiers: NSEM-140g04, 
NSEM-14b06, and NSEM-lllb04), 5 corresponding to GH18 
(identifiers: NSEM-18f06, NSEM-88dl2, NSEM-24g06, NSEM- 
96h07, and NSEM-154bl2) and 1 corresponding to GH22 (iden- 
tifier: NSEM-123b01). Only five out of nine sequences displayed 
non-truncated open read frames (OREs): NSEM-140g04 and 
NSEM-14b04 (GH16), NSEM-18f06 and NSEM-154bl2 (GH18), 
and NSEM-123b01 (GH22). Moreover, prediction of signal pep- 
tides in these sequences suggests that their putative proteins are 
secreted. According to our analysis, hypothetical complete GH 
proteins have predicted molecular masses of 40-41, 43-47, and 
15kDa and estimated isoelectric points of 6.2-7.8, 6.8-8.4, and 
4.9 (GHE16, GHE18, and GHF22, respectively). 

Sequences of interest found in the L. longipalpis EST library 
were BLASTed against NCBI-NR and identified according to 
their best hit. All analyzed sequences were identified as proteins 
belonging to families GH16, GH18, and GH22, as shown in 
Table S3. 

Amino acid sequences alignments of L. longipalpis pro- 
teins from GHE16, GHE18, and GHE22 with members from 
different insect groups exhibited highly-conserved regions (in 
Supplementary Material). L. longipalpis GH16 sequence of clone 
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NSFM-140g04 was constituted by 385 amino acid sequence and 
highly-conserved regions when compared to other insect P-1,3- 
glucanases. These conserved regions also include typical family 
GH16 catalytic glutamate residues with proton donor-acceptor 
functions and are located at positions 196 and 201. Also, the 
sequence has a putative signal peptide with cleavage sites at 
positions 19 and 20 (Figure SI). 

Putative P-glucan binding protein clones NSFM-lllbOl and 
NSFM-14b06 found in the L. longipalpis EST library were 
constituted by 268 and 370 amino acid residues, respectively. 
These sequences shared similar regions with insect P-glucan 
binding proteins, also known as Gram-negative-binding pro- 
teins (GNBPs), (Figure S2). Since NSFM-lllbOl was a trun- 
cated sequence at the 5' region, it was not possible to evaluate 
the presence of a signal peptide. On the other hand, NSFM- 
14b06 did display a signal peptide with cleavage sites between 
positions 17 and 18. Interestingly, NSFM-lllbOl and NSFM- 
14b06 lacked the catalytic glutamate residues responsible for 
hydrolysis of P-l,3-glucan otherwise present in P-l,3-glucanases 
(Figure S2). 

Alignment of family GHlS-like sequences showed that clones 
NSFM-154bl2, NSFM-88dl2, and NSFM-18f06 had a high level 
of conservation with other insect chitinases. Results of chitinases 
and chitinases-like domain organization analysis are summarized 
in Figure S3. 

Clone NSFM-96h07 consisted of a small 113 amino acid region 
located in the protein catalytic domain. It was not possible to ver- 
ify the presence of a signal peptide because 3' and 5' ORFs regions 
were truncated. 

Insect chitinases share four different conserved regions (CRl 
to CR4). Clone NSFM-154bl2 consisted of 392 amino acids and 
showed a putative signal peptide with cleavage sites between 
amino acids 19 and 20. This sequence shared similarity in CR2 
and CR3 (Figure S4). However, CRl and CR4 were not found in 
this sequence. Clone NSFM-88dl2 consisted of 296 amino acids 
and was truncated at region 5', lacking a signal peptide and CRl. 
On the other hand, CR2, CR3, and CR4 were present and CR2 
displayed the catalytic glutamate (E40) (Figure S5). 

NSFM-24g06 consisted of 201 amino acid residues. The align- 
ment of this sequence with other insect chitinases allowed identi- 
fication of CR3 and also suggests that NSFM-24g06 is incomplete 
(Figure S6). The presence of CRl and CR2 and the presence of 
catalytic residues could not be evaluated. 

Clone NSFM-96h07 showed low conservation in amino acid 
sequence in CRl and CR2 and no evidence of CR3 and CR4. 
Furthermore, NSFM-96h07 displayed a substitution of glutamate 
(E) with glutamine (Q75) at the hypothetical catalytic residue 
(Figure S7). 

NSFM-18f06 consisted of 441 residues. This sequence showed 
a signal peptide with cleavage sites between amino acids 21 and 
22. CRl and CR2 were identified within this sequence. A ser- 
ine residue (S152) was identified in CR2, unlike the all other 
sequence analyzed which had a glycine and alanine residue at 
that position. CR3 and CR4 were not identified in the alignment 
(Figure SB). 

Lysozyme-like clone NSFM-123b01 consisted of 165 amino 
acid residues and appeared to be complete as suggested by the 



presence of a methionine residue and a stop codon. Furthermore, 
SignallP showed the presence of signal peptide cleavage sites 
between amino acids 22 and 23, (Figure S9). 

PHYLOGENETIC ANALYSIS OF GH16, GH18, AND GH22 SEQUENCES 
FOUND IN THE L LONGIPALPIS ESTs LIBRARY 

To classify all L. longipapis genes found in families GH16, GH18, 
and GH22 and to predict a physiological role of these putative 
enzymes, phylogenetic trees were constructed using a Neighbor- 
joining algorithm. Subgroups considered for phylogenetic anal- 
ysis were P-glucanases and P-glucan binding proteins for family 
GH16, 1 to Vlll-chitinase groups of families GH18 (Zhang et al, 
201 1), and i and c-type lysozymes for family GH22 (Bachali et al, 
2002; Paskewitz et al, 2008). Analysis of GH16 was restricted 
to sequences from order Diptera while in GH18 and GH22 it 
considered several orders of insects. 

Phylogenetic analysis of family GH16 revealed that sequences 
found in the L. longipalpis ESTs library separately groups with 
each of the two major clades found in this family: one clade 
includes sequences with high bootstrap, all lacking the catalytic 
glutamates (e.g., NSFM-lllb04), and other clade which includes 
sequences with low bootstrap and catalytic residues (e.g., NSFM- 
140g04) (Figure lA). Moreover, the latter clade forms a para- 
phyletic group with a C. quinquefasciatus sequence annotated 
in GenBank as a Gram-negative bacteria binding protein, which 
harbors the glutamate catalytic residues. 

There was no evidence of any Pfam domain for NSFM-14b06. 
Furthermore, NSFM-lllbOl and NSFM-14b06 sequences lacked 
the catalytic site responsible for hydrolysis of P-l,3-glucan other- 
wise present in P-l,3-glucanases (Figure IB). Taking into account 
the presence and location of a catalytic site in the sequences 
within the clade of digestive enzymes, we named NSFM-140g04 
as LipClu, whereas NSFM-1 1 lb04 and NSFM-14b06 were named 
LlGBPl and L1GBP2, respectively. 

Phylogenetic analysis of GH18 divided 4 L. longipalpis 
chitinases into different clades. To perform this, we extracted 
the GH18 chitinase catalytic conserved domains (Pfam 00704), 
except for the NSFM-96h07 sequence which did not contain such 
domain. 

In our analysis, chitinase subgroups I, II, III, V, VI, VII, 
and VIII formed monophyletic groups, while group IV formed 
a paraphyletic group (Figure 2). According to clade division, 
L. longipalpis sequences were divided as follows: NSFM-154bl2 
into type VIII, NSFM-88dl2 into type II, NSFM-24g06 into type 
IV, and NSFM-1 8f06 into type V (Figure 2). After the chitinases 
classification from our filogenetics analysis, putative chitinases 
were named as follows: NSFM-154bl2 is LlChit2, NSFM-88dl2 
is LlChit3, NSFM-24g06 is LlChit4, NSFM- 96h07 is LlChitS, 
and NSFM-18f06 is LIIDGF. The name LlChitl was not used 
in this work to avoid confusion with the sequence described by 
Ramalho-Ortigao et al. (2005; LlChil). 

The phylogenetic tree of lysozymes from insects was divided 
in two clades: c-type lysozyme and i-type lysozymes, forming 
monophyletics groups in each clade (Figure 3). Positioning of 
the NSFM-123b01 sequence showed that this sequence belong- 
ing to the i-type lysozyme family (Figure 3), named in this paper 
as LlLysi. 
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FIGURE 1 I Cladogram of selected GH16 proteins from the order DIptera. 

(A) Phylogenetic tree was done using sequences from Armigeres subalbatus 
(accession nunnber AAT99011), Glossina morsitans morsitans (ABC25063), 
Drosophila melanogaster (AAF33851 ), Anopheles gambiae (ACN38130), 
Culex quinquefasciatus (AEQ27734), and Simulium vittatum (ACH56895). 
Actinophrys sol (Porifera, BAG32349) was used as outgroup. The sequences 



EXPRESSION OF p-1.3-GLUCANASES, CHITINASES, AND LYSOZYME 
DURING DEVELOPMENT OF L. LONGIPALPIS 

The expression of putative glycoside hydrolase-coding transcripts 
from families 16, 18, and 22 (P-glucanases/f5-glucan binding pro- 
teins, chitinases, and lysozyme, respectively) was evaluated at 
different development stages and dietary conditions of L. longi- 
palpis, namely larvae (L4), pupae (P), unfed male adults (UM), 
sugar-fed male adults (SM), unfed female adults (UF), sugar-fed 
female adults (SF), and blood-fed female adults (BF). Among 
GH16 sequences, Ll^Glu showed a significantly higher relative 
expression in the L4 larval stages (0.40 ± 0.09) when compared 
to the other stages of development (p < 0.05) (Figure 4A). In 
contrast, binding proteins were similarly expressed throughout all 
stages (LlGBPl, Figure 4B) or poorly expressed in the larval stage 
(L1GBP2, Figure 4C). 

LlGBPl showed a higher relative expression in P (0.89 ± 
0.05), SM (0.74 ± 0.01), and BF (0.7 ± 0.1). Relative expres- 
sion was significantly higher in P compared to L4 {p < 0.005), 
and UF and SF (p < 0.001). Furthermore, the relative expres- 
sion of SM was significantly different compared with UF and 
SF (p < 0.0001). Expression in blood-fed females was statistically 
difference compared to SF (p < 0.05) (Figure 4B). 

L1GBP2 showed higher relative expressions in the adult stage, 
specially in UM (3.4 ± 0.5), SM (4.7 ± 0.7), and BF (3.9 
± 0.5). However, expression in UM was only significantly dif- 
ferent from L4 (p < 0.05), SM expression values were signifi- 
cantly different to L4 (p < 0.01) and P (p < 0.05), while val- 
ues in BF were different from L4 (p < 0.01) and P (p < 0.01) 
(Figure 4C). 



NSFM-111 b04 and NSFM-140g04 retrieved from Lutzomyla longipalpis EST 
are named as LlGBPl and LIBeta-Glu, respectively. Bootstrap values were 
obtained by neighbor-joining method (software MEGA 5.05) using 5000 
replications. (B) Aligment of the p-1 ,3-glucanases active site sequences used 
for phylogenetic analysis showing the presence of the catalytic region (black 
box). 



Chitinases-lLke transcripts also showed different expression 
patterns throughout the sand fly developmental stages. LlChit2 
showed no significant differences between experimental groups 
(Figure 5A). 

LlChit3 showed higher relative expression in L4 (1.8 ± 
0.1), P (2.7 ± 0.6) and UM (2.6 ± 0.6). Values were signif- 
icantly different when compared to SM (p < 0.01), whQe val- 
ues in P were different from SM and UF, SF, and BF (p < 
0.05) (Figure 5B). LlChit4 showed high relative expression lev- 
els in L4 (2.8 ±0.1) and SM (2.6 ± 0.4). Expression levels 
in L4 were significantly different from UF, SF, and BF (p < 
0.001). The sugar-fed males showed significant differences com- 
pared to all adult female groups UF, SF, and BF (p < 0.05) 
(Figure 5C). 

LlChit5 showed higher relative expression in P (3.1 ± 0.3) 
compared to all groups (p < 0.05) except UM (Figure 5D). 
LIIDGF showed significant differences in relative expression in 
SM ( 16 ± 9) compared to groups L4, P, UF, SF and BF (p < 0.05) 
(Figure 5E). 

LlLysi showed higher relative expression values UM (5 ± 2). 
However, significant differences were only observed in L4 and BF 
(p < 0.05), as shown in Figure 6. 

TISSUE-SPECIFIC EXPRESSION OF ^-1,3-GLUCANASES, CHITINASES, 
AND LYSOZYME IN LARVAE OF L. LONGIPALPIS 

LipGlu showed a significantly higher relative expression level in 
larval guts when compared to other tissues and it appears to be 
gut-specific (p < 0.01, Figure 7A). On the other hand, expres- 
sion of LlGBPl was significantly higher in carcass (Figure 7B) and 
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FIGURE 2 I Cladogram of selected insect GH18 protein sequences and 
positioning of L. longipalpis GH18 (LIChit2, LIChitS, LIChit4, and LIIDGF) 
in chitinase subgroups. Sequences are from: Aedes aegypti (accession 
number XP_001657537; XP_001 656234; XP_001656233; XP_001656232; 
XP_001656231; XP_001 662588; XP_001650020; XP_001 663097; 
XP_001 656054; XP_001 655973; XP_001 654045; XP_001 663568; 
XP_001655071; XP_001 663099; XP_001660745; XP_001 660748), Anopheles 
gambiae (XP_315650; XP_315351 ; HQ456129; HQ456130; HQ456131; 
HQ456132; HQ456133; XP_308858; XP_316448; XP_307732; XP_001 238192; 
XP_310662; XP_316142; XP_314312; XP_319801 ; XP_001 688641 ; 
XP_316256; XP_001 237925), Drosophila melanogaster (NP_477298; 



NP_524962; NP_650314; NP_572598; NP_647768; NP_6n542; NP_6n543; 
EAA46011; NP_572361; NP_726022; NP_477258; NP_477257; NP_723967; 
NP_727374; NP_611321; NP_477081 ) and Tribolium castaneum 
(NP_001034516; NP_001073567; NP_001034524; XP_967813; 
NP_001036035; NP_001038094; NP_001038096; NP_001036067; XP_974461 ; 
XP_972802; NP_001036034; XP_973005; XP_973077; NP_001034515; 
XP_972719; XP_973161 ; XP_973119; XP_970191; NP_001 034517; 
NP_001038095; NP_001038092; NP_001038091. Bootstrap values were 
obtained by neighbor-joining method (software MEGA 5.05) using 5000 
replications. Branches marked with an asterisk correspond to L. longipalpis 
sequences. 



L1GBP2 showed higher expression levels in head and carcass when 
compared to the gut (p < 0.01) (Figure 7C). 

LlChit2 and LlChit3 showed no statistically significant differ- 
ences between any tissue-specific samples (Figures 8A,B, respec- 
tively). However, data suggest a slightly higher expression in 
carcass samples for LlChitS (Figure 8B). 

LlChit4 was more expressed in the gut tissue when compared 
to the head and carcasss of sand fly larvae (Figure 8C, p < 0.01). 
LlChitS was more expressed in carcass when compared to gut 
samples (p < 0.01, Figure 8D). LIIDGF showed significantly 
higher expression values in carcass when compared to other tis- 
sues {p < 0.01, Figure 8E). LlLysi also showed higher expression 
levels in carcass compared to the head or intestine (p < 0.01, 
Figure 9). 



DISCUSSION 

DIGESTION OF MICROORGANISMS IN DETRITIVORE INSECTS AND 
SAND FLIES 

Despite the wide variety of dietary sources used by insects 
in nature, feeding on decaying organic matter (i.e., plant 
debris and animal feces) is a recurrent evolutionary trait 
in several insect orders as Dictyoptera, Isoptera, Coleoptera, 
Diptera. Furthermore, some insects feed on organic mat- 
ter decomposers such as fungi and are called fungivo- 
rous. These groups include some social insects from order 
Hymenoptera so specialized that they cultivate fungi to feed 
their own colonies (Chapela et al, 1994). Some species 
of beetles, such as Dorcus rectus feeds on decaying wood 
and its associated fungi in nature and can survive under 
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FIGURE 3 I Cladogram of selected insect GH22 protein sequences (c and 
i-type lysozymes) and LILysi. Sequences are from: Anopheles darlingi 
(accession number AAB61345), Anopheles gamblae (AAC47326; AY659931 ; 
EF492429), Drosophlla melanogaster (AAF47445; AAF47453; CAA21317; 
AAF57939; AAF57940; AAF47452; AAF47451 ; AAF47450; AAF47449; 
AAF47448; CAA80225), Heliothls virescens (AAD00078), Hyalophora cecropia 



{AAA29189), Hyphantria cunea (AAA84747), Manduca sexta (AAB31190), 
Musca domestica (PC4062), Thatoma infestans (ABI94387), Thatoma 
brasiliensis (AAU04569), Rhodnlus prollxus (ABX11554; ABX11553), Culex 
tarsalis (ACJ64375), Aedes aegypti (CAC19819). Bootstrap values were 
obtained by neighbor-joining method (software MEGA 5.05) using 5000 
replications. 



laboratory conditions fed merely on fungi (Tanahashi et al., 
2009). 

L. longipalpis is a Dipteran from family Psychodidae and 
hematofagous females of this species are the main vector of vis- 
ceral leishmaniasis in the New World. However, little is known 
about the feeding habits of their larvae in their natural envi- 
ronment. In experimental conditions, larval L. longipalpis and 
L. intermedia developed better when fed on a diet of fungi- 
rich humus which mimics larval substrate in the wilderness 
(Wermelinger and Zanuncio, 2001). Recently, breeding sites for 
sand fly larvae were described in the forest located at Amazonas 
State (Alencar et al., 2011) and in urban and peri urban areas in 
Southeast Brazil (Casanova et al., 2013). In the forest, sand fly 
larvae seem to be associated with the humus-rich soil near tree 



bases, and in urban and peri-urban areas the preferential breeding 
sites seem to be the soil at chicken sheds. Despite the limitations 
of these studies, it is quite relevant to perceive the strict associa- 
tion of sand fly larvae with microorganism-rich decaying organic 
matter in the wild. 

Considering that, it is expected that sand fly larvae could 
exploit microorganisms as bacteria and fungi as a nutritional 
source. In this respect, their digestive enzymes must be capa- 
ble of chemically disrupt the cell walls of these microorganisms, 
in order to release intracellular molecules as nucleic acids, pro- 
teins and storage sugars and lipids, which are essential to the 
insect development and metabolism. The main digestive enzymes 
related to fungal and bacterial cell disruption in insects are P-1,3- 
glucanase, chitinase, and lysozyme. Digesting P-l,3-glucanases 
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FIGURE 4 I Relative expression of p-1,3-glucanases (LipGlu) (A) and 
p-glucan binding proteins (LIGBP1 and LIGBP2) (B,C, respectively) in 
different stages of development of Lutzomyia longipalpis as determined 
by Multiplex RT-PCR (Ribosomal protein 60 used as constitutive gene). 

The experimental groups used in the experiments were larvae (L4), pupae 



(P), unfed male adults (UM), sugar-fed male adults (SM), unfed female adults 
(UF), sugar-fed female adults (SF), and blood-fed female adults (BF). Bars 
represent the mean ± SB from 3 independent experiments. The statistics 
were done using ANOVA with Turkey post-test and student's f-test. Asterisks 
indicate statistically significant differences in the samples. 



from GHF16 were already described in detritivore insects from 
orders Dictyoptera (Genta et al., 2003), Coleoptera (Genta et al, 
2009), Isoptera (Lucena et al., 2011). Digestive chitinases from 
GHF18 were described in detritivore Coleoptera (Genta et al, 
2006), and digestive Lysozymes (GHF22) are a common feature 
of Dipteran Cyclorrapha (Terra and Ferreira, 2005). All these 
enzymes were found in a recent report on sand fly larval gut activ- 
ities (Moraes et al., 2012) but, contrarily to the examples above, 
sand fly enzymes could not be molecularly identified by tradi- 
tional biochemical techniques, due to the minimal size of these 
insects. In spite of that, our bioinfomatic and expression analysis 
of selected transcripts in L. longipalpis suggest that these larvae 
use GHF16 and GHF18 in their digestion, as observed in other 
insect orders. 

FUNCTIONAL SPECIALIZATION OF GHF16 MEMBERS IN L 
LONGIPALPIS 

The in silico search for GH-like sequences in a L. longipalpis 
EST library showed the presence of three clones highly simi- 
lar to family GH16 proteins (NSFM-14g04, NSFM-lllb04, and 
NSFM-14b06). Previous sequence comparisons and phylogenetic 
analysis supported an evolutionary relationship between P-1,3- 
glucanases and P-glucan binding proteins from GHF16 (Pauchet 
et al, 2009; Bragatto et al, 2010; Hughes, 2012). The most evi- 
dent features which distinguish these two functional groups are 



(1) the presence of catalytic glutamates in a conserved region 
which correspond to the active site in P-l,3-glucanases and (2) the 
presence of a conserved 100 amino acid N-terminal extension in 
P-glucan binding proteins. It has been proposed that the animal 
P-l,3-glucanase ancestral gene suffered a duplication before the 
differentiation of arthropods and molluscs (Bragatto et al, 2010) 
and, in this respect, insects should bear at least two copies of genes 
from GHF16. 

Pauchet et al. (2009) also divided clades of Lepidopteran 
GHF16 sequences according to the absences or presences of 
catalytic residues. The authors suggested that such division hap- 
pened through a duplication event of a gene in a common ances- 
tor. This originated two major groups with different functions 
in insects: digestion and immune signaling. Recently, Hughes 
(2012) showed that GHF16 Pathogen Recognition Receptors 
(PRRs, which includes {5-1,3-glucan binding proteins) and P-1,3- 
glucanase sequences are present in orders Coleoptera, Diptera 
and Lepidoptera, and suggested that this gene duplication event 
occurred before divergence of these holometabolous orders. The 
author also exhibited that only PRR-like sequences are found 
within Exopterigota (hemimetabolous), while only glucanase- 
like sequences are found in Isoptera. Such evidence supports 
the hypothesis that these two subfamilies underwent events of 
gene duplication before the origin of holometabolous insects (i.e., 
during the Carboniferous, approximately 300 million years ago). 
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FIGURE 5 I Relative expression of GHF18 members in different stages of 
development of Lutzomyia longipalpis as determined by Multiplex 
RT-PCR (Ribosomal protein 60 used as constitutive gene). Chitinases 

amplified were: (A) LIChit2; (B) LIChitS; (C) LIChit4; (D) LIChitB, and (E) 
LIIDFG. The experimental groups used in the experiments were larvae {L4), 



pupae (P), unfed male adults (UM), sugar-fed male adults (SM), unfed female 
adults (UF), sugar-fed female adults (SF), and blood-fed female adults (BF). 
Bars represent the mean ± SE from 3 independent experiments. The 
statistics were done using ANOVA with Turkey post-test and student's f-test. 
Asterisks indicate statistically significant differences in the samples. 



Based on the criteria described above, NSFM-14g04 was clas- 
sified as a P-l,3-glucanase (LipClu) while NSFM-lllb04 and 
NSFM-14b06 were assigned to the f5-glucan binding protein 
group (as LlBGPl and L1BGP2). Alignment of these sequences 
with other insect GH16 showed that clone NSFM-14g04 con- 
tains the two glutamate residues which are important for catalysis. 
In contrast, NSFM-11 lb04 and NSFM-14b06 lack these residues, 
which suggests that these putative proteins do not have enzy- 
matic activity. The presence of the N-terminal region typical for 
P-glucan binding proteins was not confirmed in LlGBPl, because 
this sequence is truncated at the N-terminal, but was confirmed 
in L1GBP2. LipGlu, as expected, does not contain this extension, 
which was confirmed as LlfiGlu appears as a complete coding 
sequence with a putative signal peptide. The presence of a puta- 
tive signal peptide in LipGlu and L1GBP2 sequences is coherent 



with the observations that some digestive insect P-l,3-glucanases 
foUow the exocytic route for secretion (Bragatto et al., 2010) 
and that P-glucan binding proteins are soluble proteins secreted 
to the hemolymph, where they interact with members of the 
prophenoloxidase activating cascade (Lee et al., 2004). 

The expression pattern of GHF16 members in L. longipalpis 
corroborated the classification based on sequence features. LlfiGlu 
was more expressed in larvae and, at this stage, its expression 
is gut-specific. Gut-specific P-l,3-glucanases from GHF16 were 
already described in Lepidoptera (Pauchet et al, 2009; Bragatto 
et al., 2010), Coleoptera (Genta et al, 2009), and Isoptera (Bulmer 
et al., 2009). To our knowledge, this is the first description of 
this kind of protein sequences in Diptera. The role of insect 
gut P-l,3-glucanases is still controversial, as they were impli- 
cated in digestion of fungi and plant cell wall polysaccharides in 
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some insects (Dictyoptera, Genta et al., 2003; Orthoptera, Genta 
et al., 2007, Coleoptera, Genta et al., 2009) but in recognition 
of gut-pathogens in others (Lepidoptera, Pauchet et al., 2009; 
Isoptera, Bulmer et al, 2009). Considering the detritivore habit 



of L. longipalpis larvae in the laboratories, with the ingestion of 
significant amounts of fungal tissue in the food (Moraes et al., 
2012), we are tempted to suppose that the main role of Ll^Glu 
is the disruption of ingested fungal cells, but more functional 
studies are necessary to confirm this hypothesis. 

At the same time, LlGBPl and L1GBP2 showed low levels of 
expression in the larvae and, at this stage, low levels of expression 
in the gut. In this respect, a role for these genes in larval digestion 
can be ruled out. Besides that, the expression patterns of LlGBPl 
and L1BGP2 strongly suggest that they are involved in defense 
against pathogens, as they are expressed in tissues and stages more 
susceptible to infection. Both genes are highly expressed in males 
fed with sugar and blood fed females. The ingestion of bacterial 
pathogens by sand flies in sugar meals is a current topic of investi- 
gation (Telleria et al., 2013), as well as the multiplication of these 
micro organisms inside the blood meal in the sand fly gut (Diaz- 
Albiter et al., 2012). The fact that sand flies mount strong defenses 
in the gut during these physiological conditions suggests that they 
are fighting pathogens and, in fact, interfere with this phenomena 
result in activation of the prophenoloxidase cascade and insect 
death (Diaz-Albiter et al, 201 1). 

Nevertheless, there are subtle differences in the expression pat- 
terns between LlGBPl and L1GBP2. LlGBPl is highly expressed 
during the pupal stage, as well as L1GBP2 is highly expressed 
in unfed males. In this respect, these genes behave like their 
putative homologs in Drosophila, where GNBPl (CG6895-PA), 
GNBP2 (CG4144), and GNBP3 (CG5008) are not expressed in 




Development Stages 



FIGURE 6 I Relative expression of lysozyme (LILysi) in different stages of 
development of Lutzomyia longipalpis as determined by Multiplex 
RT-PCR (Ribosomal protein 60 used as constitutive gene). The 

experimental groups used in the experiments were larvae (L4), pupae (P), 
unfed male adults (UM), sugar-fed male adults (SM), unfed female adults (UF), 
sugar-fed female adults (SF), and blood-fed female adults (BF). Bars represent 
the mean ± SB from 3 independent experiments. The statistics were done 
using ANOVA with Turkey post-test and student's Most. Asterisks indicate 
statistically significant differences in the samples. 




FIGURE 7 I Relative expression of GHF16 (A) p-1,3-glucanase (LipGlu) 
and (B) p-glucan binding proteins LIGBP1 and (C) LIGBP2 in 
different tissues of Lutzomyia longipalpis L4 larvae as determined 
by Multiplex RT-PCR (Ribosomal protein 60 used as constitutive 



gene). Bars represent the mean ± SB from 5 independent 
experiments. The statistics were done using ANOVA with Turkey 
post-test and student's Mest. Asterisks indicate statistically significant 
differences in the samples. 
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FIGURE 8 I Relative expression of GHF18 chitinases and chitinase-like 
proteins in different tissues of Lutzomyia longipalpis L4 larvae as 
determined by Multiplex RT-PCR (Ribosomal protein 60 used as 
constitutive gene). Chitinases amplified were: (A) LIChit2; (B) LIChitS; (C) 



LIChit4; (D) LIChitB and (E) LIIDGF Bars represent the mean ± SE from 5 
independent experiments. The statistics were done using ANOVA with 
Turkey post-test and student's f-test. Asterisks indicate statistically significant 
differences in the samples. 



the midgut at any stage, with higher expression levels in the 
embryo (GNBP2 and 3), pupae (all GNBPs) and adults (GNBP2 
and 3), and preferential expression in the carcass or head (St. 
Pierre et al., 2014). It is possible that LlGBPl is related to the pre- 
vention of infections during the pupal stage, which is extremely 
sensitive to infections. Besides that, it has been shown that unfed 
sand flies carry a significant amount of bacteria from the lar- 
val stage (Sant'Anna et al, 2014), which suggests that L1GBP2 
could be involved in the specific protection against these recal- 
citrant microorganisms. The trans stadial passage of bacteria is a 
well-documented phenomenon in sand flies (Volf et al., 2002), 
and may have important implications in the development of 
strategies for the blocking of vectorial transmission of pathogens 
based in paratransgenesis. However, more functional studies are 



needed to confirm the roles of LlGBPl and L1GBP2, especially 
after challenge with insect pathogens. 

FUNCTIONAL SPECIALIZATION OF GHF18 MEMBERS IN L 
LONGIPALPIS 

Five sequences from GHF18 were retrieved from the L. longi- 
palpis EST library (NSFM-96h07, NSFM-154bl2, NSFM-88dl2, 
NSFM-24g06, and NSFM-18f06). According to amino acid sim- 
ilarities, domain compositions, and phylogenetic analysis, insect 
proteins from GHF18 are classified into 8 groups, I- VIII (Zhang 
et al, 2011). A key feature in these proteins is the presence of 
a catalytic glutamate in the sequence DWEYP at the consensus 
region 2 (CR2), which is used as a marker for enzymatic activ- 
ity. GHF18 proteins without this residue are named chitinase-like 
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FIGURE 9 I Relative expression of i-type lysozyme (LILysi) in different 
tissues of Lutzomyia longipalpis L4 larvae as determined by IVIultiplex 
RT-PCR (Ribosomal protein 60 used as constitutive gene). Bars 
represent the mean ± SE from 5 independent experiments. The statistics 
were done using ANOVA with Turkey post-test and student's f-test. 
Asterisks indicate statistically significant differences in the samples. 



proteins and considered devoid of hydrolytic activity. Besides 
that, functional studies were able to incriminate groups I and II 
of insect chitinases in chitin hydrolysis during molting, group 
III in the distension of wings and abdomen during morpho- 
genetic development, and group IV as digestive enzymes. Group 
V is devoid of catalytic activity and is referred as Imaginal Disk 
Growth Factors, promoting cell proliferation in imaginal disks. 
Groups VI -VIII were described based only in sequence attributes, 
with no clear functional role assessed to this moment. In gen- 
eral, chitinases from all groups but IV are expressed in all stages 
of development, whUe the expression of group IV chitinases is 
restricted to specific stages. Many group IV chitinases are larval- 
specific, and gut-specific genes, being considered primarily as 
digestive enzymes. Some are involved in digestion or turnover of 
the peritrophic membrane (Ramalho-Ortigao and Traub-Cseko, 
2003), but others, which lack the characteristic C-terminal Chitin 
Binding Domain, are believed to act in disruption of fungal cell 
walls without affecting the PM structure (Genta et al., 2006). 

The domain organization of all L. longipalpis GHF18 
sequences present in the ESTs database could not be assessed, 
because only two of them (NSFM-154bl2 andNSFM-18f06) code 
for complete ORFs. These complete sequences contain a putative 
N-terminal signal sequence, which is a common feature in GHF18 
insect chitinases (Zhang et al., 2011). This is coherent with chiti- 
nase putative secretion to the molting fluid, intestinal lumen, or 
to the extracellular matrix where it exerts its action on imaginal 
disk cells (Arakane and Muthukrishnan, 2010). 

In spite of lacking complete N- or C-terminal ends, all GHF18 
proteins from L. longipalpis were assigned to functional groups of 
insect chitinases II, IV, V, or VIII, being named LlChit2 (group 
VIII), 3 (group II), 4 and 5 (Group IV) and LIIDGF (group 
V). From these sequences, only LlChit2 and LlChit3 contain the 
conserved catalytic glutamate at CR2, being probably true chiti- 
nases. LIIDGF is a chitinase-like protein, and it was not possible 
to check for the presence of the catalytic glutamate in the case 
of LlChit4, because CR2 was missing in this sequence. In spite 



of that, LlChit4 is probably an active enzyme, as it groups with 
other true chitinases from group IV and contain the conserved 
PFAM domain Glyco_hydro_18 (PF00704). LlChitS also ahgned 
with proteins of group IV, but was not included in the phyloge- 
netic analysis because the PFAM domain PF00704 was missing in 
its sequence. Because of that, the presence of the catalytic gluta- 
mate in region CR2 could not be really assessed, as its apparent 
substitution in LlChitS sequence could be the result of the forced 
alignment of its truncated sequence. Another possibility is that 
LlChit5 use an alternative basic hydrolysis mechanism, based only 
in the nucleophUic catalysis from the acetamido group of the sub- 
strate (Hashimoto et al., 2000). More evidence is necessary to 
ascertain the catalytic nature of LlChit4 and LlChitS. 

Chitinases from subgroups II, V, and VIII are typically 
expressed in all stages of insect development (Zhang et al., 201 1). 
The expression pattern observed for LlChit2, 3, 5, and LIIDGF 
is consistent with this behavior. Indeed, LlChit2, 3, and 5 reach 
their highest expression levels in the pupal stage, and LIIDGF is 
more expressed in adult males. Besides that, in the larvae these 
genes are more expressed in the carcass (LlChit3, 5, and LIIDGF) 
or equally expressed in all larval tissues (LlChit2). These data 
strongly suggest that LlChit2, 3, 5, and LIIDGF are not involved 
in the digestion of sand fly larvae. 

Interestingly, LlChit4, which is related to typical insect diges- 
tive chitinases from group IV, is more expressed in larvae and 
sugar fed males. Type IV chitinases contain a signal peptide, a 
single catalytic domain and mostly no CBDs (Genta et al., 2006). 
Additionally, LlChit4 is more expressed in the gut of larvae, which 
strongly suggest that this transcript correspond to the larval diges- 
tive chitinase of L. longipalpis. The higher expression in sugar 
fed males could be related to the proliferation of microorgan- 
isms in the gut of these insects, but this hypothesis still needs 
confirmation. 

Similar results were observed by Khajuria et al. (2010), who 
demonstrated by phylogenetic analysis that a chitinase (OnCht) 
predominantly expressed in the gut of larval Ostrinia nubilalis was 
a type IV chitinase. A digestive chitinase was also found in the 
midgut of T. molitor larvae (Genta et al., 2006, TmChi). 

In several hematophagous insects, digestive chitinases have an 
extremely important role in the degradation of type I peritrophic 
matrix (PM I). Zhang et al. (201 1 ) showed expression of numer- 
ous chitinases in Anopheles gambiae, among them, a chitinase 
mainly expressed in the gut of adults (AgChtS). 

Previous studies made by Ramalho-Ortigao and Traub-Cseko 
(2003) demonstrated the expression of chitinase (LlChil) in the 
midgut of adult female L. longipalpis, whose expression seems to 
be induced after blood feeding. However, so far the expression 
of chitinases in larvae of L. longipalpis has not been demon- 
strated and studies to pursue a chitinase important for digestion 
of microorganisms present in the diet of the larvae are still scarce. 
The data obtained suggest that LlChit4 probably has that role. 

To our knowledge, this is the first description of a larval gut 
chitinase sequence in sand flies. Larval gut chitinase activities 
were already described in L. longipalpis (Moraes et al., 2012) 
and Aedes aegypti (Souza-Neto et al., 2003), these enzymes being 
probably involved in the digestion of chitin rich structures as 
the fungal cell wall. Gut-specific expression of chitinase genes in 
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Diptera was already described in Anopheles gambiae larvae (Zhang 
et al, 2011), where group IV chitinases AgChtS and AgChtl3 are 
majorly expressed only in the midgut at this stage. It is interesting 
to notice that in some cases dipteran larval digestive chitinases 
seem to be coded by genes different from their adult counterparts, 
a pattern already observed for other gut genes, as trypsins and 
peritrophins (Venancio et al., 2009). 

ANALYSIS OF THE GHF22 SEQUENCE FOUND IN ADULT ESTs LIBRARY 
OF LUTZOMYIA LONGIPALPIS 

Lysozymes are present in a wide variety of organisms (from 
viruses to plants and animals) and can be classified based on 
sequence similarity into 5 groups: type c, type g, type i, plant, 
and bacterial lysozymes (Bachali et al., 2002). Insect lysozymes 
are commonly members of family GH22, being involved in the 
arthropod innate immune response against bacteria but also in 
the digestion of food microorganisms in detritivores (JoUes and 
JoUes, 1984; Hultmark, 1996). The most well characterized insect 
digestive lysozyme is the gut activity from Musca dotnestica lar- 
vae. This enzyme is a c-type lysozyme (Can^ado et al., 2007) 
with specializations in its structure which are related to its opti- 
mum pH at the midgut strong acidic conditions (Can^ado et al., 
2010), which is a typical feature of Diptera Cyclorrapha (Terra 
and Ferreira, 1994). As the midgut luminal pH of sand fly larvae 
is strongly alkaline (Do Vale et al, 2007) it could be possi- 
ble that their digestion of bacteria rely on a different type of 
lysozyme. 

In recent years studies have revealed the existence in insects of 
a second family of lysozymes, type i, for which at least 20 exam- 
ples have been identified in the phylum Mollusca, Nematoda, 
Annelida, and Arthropoda (Ito et al, 1999; Zavalova et al, 2000; 
Bachali et al, 2002; Paskewitz et al, 2008). 

Therefore, considering the possible role of lysozyme in the 
degradation of dietary microorganisms in sand flies, a search 
was also made for GHF22 members in the EST database of 
L. longipalpis. Our analysis revealed only one transcript similar to 
insect lysozymes, called LlLysi (NSFM-123b01), whose sequence 
was complete including a putative N-terminal signal peptide. 
Moreover, comparison of this sequence in the Pfam database 
showed that LlLysi contains a destabilase domain. 

Similarity to Destabilase is common among i-type lysozymes. 
The leech destabilase has such a name from their ability to break 
or destabilize fibrin blood clots in the ingested blood (Zavalova 
et al., 2000). Using this relationship, the leech destabilase was used 
as a pattern to pursue i-lysozymes in the genome of Anopheles 
gambiae (Paskewitz et al., 2008). 

Aiming to classify the lysozyme found in the EST library 
of L. longipalpis, we performed a phylogenetic analysis using 
the neighbor-joining algorithm. In this analysis it was possi- 
ble to demonstrate the separation of c- or i-lysozymes in two 
monophyletic groups. LlLysi was classified in this analysis as i- 
lysozyme. To our knowledge, this is the first description of this 
type of lysozyme in sand flies. 

LlLysi was expressed throughout all developmental stages of 
L. longipalpis, expression being the highest in unfed males. This 
could be related to the presence of high amounts of bacteria in 
recently emerged adults of L. longipalpis (Sant'Anna et al., 2014). 



Besides that, a very low level of LlLysi expression was found in 
L4 larvae, being concentrated in the carcass, excluding the pos- 
sibility of LlLysi being a larval digestive enzyme. Similar studies 
were done in the mosquito A. gambiae, which produces two i- 
type lysozymes, called Lys i-1 and i-2. Lys i-1 was expressed in 
all stages of mosquito development, being more expressed in the 
ovaries of females, followed by Malpighian tubules and fat body. 
Lys i-2 was more abundant in the fat body of adults (Paskewitz 
et al, 2008). 

Thus, in sand flies, the lysozyme gene involved in the larval 
digestion of bacteria, which corresponds to the activity described 
by Moraes et al. (2012) is different from the i-lysozyme transcripts 
represented in adult EST databases. Our studies raise the interest- 
ing possibility that sand fly larval gut lysozyme could be in fact 
more related to the c-lysozymes from other Diptera, but more 
studies are necessary aiming the molecular identification of this 
enzyme. 

CONCLUSION 

The comparative analysis of sequences present in EST databases as 
well as the determination of expression patterns during develop- 
ment and the tissue location of transcripts from families GH16 
and GH18, allowed us to identify the genes related to the P- 
1,3-glucanase and chitinase which are involved in digestion of 
microorganisms in larvae of L. longipalpis. However, the sequence 
of lysozyme present in adult EST databases (LlLys, family GH22) 
is apparently a type i Lysozyme involved in immunity, having no 
relationship to larval digestion of bacteria. To our knowledge, 
this is the first molecular identification of digestive enzymes in 
the larvae of sand flies, which are important vectors of several 
diseases, including Leishmaniasis. These findings may have a pro- 
found impact on the understanding of the biology of these disease 
vectors, as well as for the development of new strategies for vector 
control. 
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